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Crossover from commensurate to incommensurate antiferromagnetism in 
stoichiometric NaFeAs revealed by single-crystal 23 Na, 75 As-NMR experiments 
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We report results of 23 Na and 75 As nuclear magnetic resonance (NMR) experiments on a 
self-flux grown high-quality single crystal of stoichiometric NaFeAs. The NMR spectra revealed 
a tetragonal to twinned-orthorhombic structural phase transition at To = 57 K and an anti- 
ferromagnetic (AF) transition at Taf = 45 K. The divergent behavior of nuclear relaxation 
rate near Taf shows significant anisotropy, indicating that the critical slowing down of stripe- 
type AF fluctuations are strongly anisotropic in spin space. The NMR spectra at low enough 
temperatures consist of sharp peaks showing a commensurate stripe AF order with a small 
moment ~0.3 fiB. However, the spectra just below Taf exhibits highly asymmetric broadening 
pointing to an incommensurate modulation. The commensurate-incommensurate crossover in 
NaFeAs shows a certain similarity to the behavior of SrFe2As2 under high pressure. 
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Stripe-type antiferromagnetic (AF) order has been 
commonly observed on the undoped or low-pressure re- 
gion in the phase diagrams of iron-based superconduc- 
tivities. Thus AF fluctuations are often considered the 
driving force for superconducting (SC) paring. 1 ' For ex- 
ample, the s±-wave pairing state, where the order param- 
eter changes sign between the nested electron and hole 
Fermi surfaces, can be stabilized by the interband AF 
fluctuations. 2 ' Since the geometry of the Fermi surfaces 
depends on materials and changes with doping, the nest- 
ing wave vectors may deviate from the commensurate 
(C) conditions in some cases, leading to incommensu- 
rate (IC)-AF order, or in other words, spin-density-wave 
(SDW) order in a narrow sense. 

However, all the AF orders of the undoped sto- 
ichiometric compounds so far are known to be C- 
type. 3 ' IC order has been reported for some non- 
stoichiometric or chemically doped materials. Neutron 
scattering experiments on the lightly electron-doped sys- 
tems Fei+aTei-^Sez reported short-ranged IC order 
characterized by broad magnetic peaks off-centered from 
the commensurate wave vectors. 4 ' 5 ' IC-SDW order has 
been also reported for Ba(Fe 1 _ a; Co a ;)2As2 by NMR 6 ' and 
Mossbauer 7 ' experiments. However, microscopic disorder 
in these materials should bring distribution in the mag- 
nitude of ordered moments, which might be difficult to 
distinguish from incommensuration due to Fermi surface 
nesting. Therefore, finding a stoichiometric material ex- 
hibiting IC-SDW would mark important advance in our 
understanding of magnetism in iron-pnictide materials. 

Among different types of iron-based superconductors, 
122-type AFe2As2 (A=Ca, Sr, Ba) have been attracting 
strong interest because large single crystals with stoi- 
chiometric composition can be grown. Recently Ill-type 
materials, LiFeAs and NaFeAs (Fig. la), have also been 



found suitable for crystal growth with alkaline-rich self 
flux. 8 ' 9 ' LiFeAs is a bulk superconductor with the tran- 
sition temperature T$c of 17 K, 10 ' while NaFeAs shows 
only filamentary superconductivity. 11 ' 12 ' NaFeAs shows 
successive phase transitions: tetragonal-to-orthorhombic 
structural transition at To ~ 50 K and AF transition 
at Taf ~ 40 K. 8,12 ' The transition temperatures, which 
are significantly lower than those for 1111 and 122 ma- 
terials, 3 ' seem to depend sensitively on self-doping due 
to Na deficiency or synthetic methods. 12 ' In this letter, 
we report observation of IC-AF order and a temperature- 
driven IC-C crossover in stoichiometric NaFeAs by NMR 
experiments on 23 Na and 75 As nuclei. 

The crystals of NaFeAs were grown by the self-flux 
method with starting elemental ratio of 2:1:2. The mate- 
rials were put into alumina crucible, sealed in a doubled- 
wall quartz tube in argon atmosphere, and heated to 
900° C. Crystals grew during subsequent cooling down 
to 400°C in three days, then to room temperature in one 
day. We confirmed 1:1:1 atomic content within 2% by 
energy dispersive X-ray spectroscopy. Resistivity, mag- 
netization and specific heat measurements detected the 
successive transitions at To = 57 K and Taf = 45 K. 
The values of To and Taf of our crystal, also confirmed 
by the NMR experiments as described below, are sig- 
nificantly higher than the previous reports, 11,12 ' indi- 
cating less self-doping. The c-axis length measured by 
single-crystal x-ray diffraction was 7.05 A at room tem- 
perature, which is longer than the previous reports. The 
crystal (5 x 3 x 0.3 mm 3 ) was handled in purified ar- 
gon and covered with grease before exposed to air for 
mounting on the NMR probe equipped with a double- 
axis goniometer. The NMR spectra were acquired us- 
ing Fourier step summing technique. The Fourier trans- 
formed spectrum of spin-echo was accumulated at shifted 
frequencies while the frequency or magnetic field was 
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Fig. 1. (Color online) (a) Crystal structure and stripe- type anti- 
ferromagnetic order (red arrows) of NaFeAs. The arrows at As 
and Na nuclei indicate the hyperfine field in the antiferromag- 
netic state, (b) Field-swept 75 As-NMR spectra at 48.31 MHz. 
Below To, the satellite lines for H X c splits into two sets due 
to twining. Below Tj^p, the center line for H \\ c splits sym- 
metrically by the hyperfine fields parallel to the c axis at the 
As sites. The dotted line shows the unshifted resonance field. 
The full widths at the half maxima (enclosed by the arrows) are 
narrower than those for 122 parent compounds, 13, 14 ) indicating 
good stoichiometry of the sample, (c) Temperature dependence 
of the intensity at the peak position of the paramagnetic spec- 
trum multiplied by T for H \\ c . No hysteresis was observed 
across Taf 5 except for the small instrumetal delays. 



being swept. NMR relaxation rates were deter- 

mined by the inversion recovery method. Good fitting 
to the theoretical recovery curve for spin 1 — 3/2 nuclei, 
{exp(-t/Ji) + 9exp(-6t/Ti)}/10, was obtained in the 
whole temperature range. 

The NMR spectra are shown in Fig. 1(b) for 75 As nu- 
clei and in Fig. 2 for 23 Na nuclei. Both As and Na occupy 
the same type of site, the 2c site in the high-T tetrag- 
onal PA/nmm structure and the Ag site in the low-T 
orthorhombic Cmma structure. Since both nulcei have 
spin 3/2, the NMR spectra consist of quadrupole-split 
three lines with the resonance frequencies (here up to the 
first-order perturbation of the quadrupole interaction) . 



/res = /Jo7N#cff + iSv {% = -1,0,1), 



(1) 



where 7n is the nuclear gyromagnetic ratio (2ir x 
7.29019 MHz/T for 75 As and 2tt x 11.26226 MHz/T for 
23 Na), H e g = H cxt + Hhf is the sum of the external field 
and the magnetic hyperfine field, and 
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Fig. 2. (Color online) Frequency-swept 23 Na-NMR spectra at 
4.25 T (a) for H || c at different temperatures and (b) for H _L c 
along various directions in the ao-plane at T = 50 K. The spectra 
for H || c show incommensurate modulation immediately below 
Taf ■ The satellite lines are doubled for H _L c due to twinning in 
the orthorhombic structure. Their angular dependence in the aft- 
plane is well reproduced by eq. (2) with r\ = 0.027 (dashed lines), 
(c) A color plot of the satellite intensity as a function of tem- 
perature for H || a or 6, showing the tetragonal-to-orthorhombic 
structural transition at Tq = 57 K. 



is the quadrupole splitting. Here, (9, <ft) specifies the po- 
lar angle of if e ff with respect to the c-axes and 77 = 
\v aa — v hh \j\v cc \ is the asymmetry parameter of the elec- 
tric field gradient (EFG) tensor v aa . In the tetragonal 
PA/nmm structure, 77 must be zero. While the center 
line (i = 0) is determined only by the magnetic hyper- 
fine field, therefore, provides direct information on the 
spin structure, the satellite lines (i = ±1) are affected by 
EFG which is a sensitive probe for structural transitions. 

The 75 As-NMR spectrum at 70 K (Fig. lb) shows three 
lines for both H \\ c and H _L c as expected. At 50 K 
below Tq the satellite lines are doubled for H \\ a or 
H || 6 in the orthorhombic notation. This has been ob- 
served also in 122-type compounds and ascribed to non- 
zero 77 and twinned domains in the orthorhombic struc- 
ture. 13 ' 14 - ) Details of the transition were investigated by 
23 Na NMR (Fig. 2) providing better S/N ratio. The an- 
gular dependence of the satellite resonance frequencies 
in the a6-plane at 50 K (Fig. 2b) is well reproduced by 
Eq. (2) with 77 = 0.027 (dashed lines). The onset tem- 
perature for the splitting of the satellite lines (Fig. 2c) 
agrees with Tq = 57 K determined by the resistivity 
measurement. 

The AF transition is marked by vanishing of the NMR 
line for H \\ c near the paramagnetic resonance field 
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below Taf = 45 K (Fig. lc for 75 As). Absence of hys- 
teresis indicates the transition is continuous. At suffi- 
ciently low temperatures, each of three lines split sym- 
metrically into two lines for H || c, but no line splitting 
occurs for H _!_ c [see the spectra at T = 20 K for 75 As 
(Fig. lb) and 23 Na (Fig. 2a)]. This indicates staggered 
hyperfine fields at the As/Na sites along the c-direction 
H^f = (0, 0, ±A) with very small distribution of A. Fol- 
lowing the analysis on the similar results for 122 com- 
pounds, 13, 14 ' we conclude that the hyperfine field is gen- 
erated by a stripe AF order with the q- vector (1,0, ^)o 
(defined on the orthorhombic Cmna unit cell) and the 
AF moment along the a-axis (Fig. la). The off-diagonal 
element (B ac ) of the hyperfine coupling tensor relates the 
AF moment a and A as A = 4_B ac <7. 13 ' Assuming the 
same value B ac — 0.4 ~ 0.5 T//ie obtained for 122 com- 
pounds, 13 ' 14 ' the AF moment in NaFeAs is estimated to 
be 0.28 ~ 0.34 /j,b at 20 K. This is close to the value re- 
ported by Yu et al. also from NMR experiments 11 ' but 
significantly larger than the value 0.09±0.04 fi^ obtained 
by neutron experiments. 12 ' 

However, the spectral shape exhibits unusual temper- 
ature dependence near Taf- The narrow six lines of 
the 75 As NMR spectrum at 20 K develop asymmetric 
distribution at 40 K with a sharp edge at the maxi- 
mum value of A (Fig. lb). The sharp upper cutoff is 
unlikely to be caused by extrinsic disorder or inhomo- 
geneity. This spectral broadening is more obvious in the 
23 Na spectra immediately below Taf (Fig. 2a), indicat- 
ing incommensurate modulation of the hyperfine field. 
However, the spectral shape changes gradually to rather 
normal discrete peaks at lower temperature, pointing to 
an incommensurate-commensurate crossover. Since the 
broadening is observed only for H || c, the hyperfine 
fields at the As/Na sites should still be almost parallel 
to the c axis, Hm = (0,0, ±A). The spectral shape then 
represents the spatial distribution of A or the AF mo- 
ments. Therefore, the NMR spectra near Taf indicate a 
longitudinal SDW order at some incommensurate wave 
vector. 

To illustrate the real-space image of moment distribu- 
tion, we assume a periodic modulation of A between zero 
and A max along x over the wave length A. The spectral 
intensity for H \\ c as a function of A is given as 



J(±A) oc 



dA(x) 



dx 



(3) 



Integrating the experimental spectra, therefore, gives the 
relation between A and x, 



A ,f A I(A')rfA' 

2/ A_ /(A , )dA ,- 



(4) 



By converting A to the AF moment, the spatial distri- 
bution of AF moments reproduced by 23 Na-NMR is ob- 
tained as shown in Fig. 3 at various temperatures. At 
low temperatures, the AF moment is nearly uniform over 
most of the space and the rapid change occurs only in a 
narrow region. Such a profile can be regarded as usual 
C-AF domains separated by domain walls. Near Taf, 
however, the domain wall gets significantly broadened, 
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Fig. 3. (Color online) Spatial distribution of Fe moments in the 
AF state of NaFeAs, assuming the simple periodicity (see text). 
The panel (a) corresponds to the region 0<©<7r/2in the 
schematic in (b). (b) The temperature dependence of the SDW 
patterns, illustrated by normalizing the data in (a). 



making the profile closer to a sinusoidal SDW. The mi- 
croscopic mechanism for such IC-C crossover is still an 
open question. 

The IC-SDW order at (1 + 5,0, \)o or (M, §)o breaks 
cancellation of the a-component of the hyperfine field, 13 ' 
therefore, should broaden the NMR spectra for H _L 
c. The full-width of the spectrum for H ± c due to 
incommesuration is given by ~ 8B aa atan(7r<5), where 
B aa (~ 75 A ab /4) is the diagonal element of the hyper- 
fine coupling tensor determined from the Knight shifts 
described below. The observed width (0.02 T at 43 K) 
then puts an upper limit for <5, \5\ < 0.004. Such tiny 
deviation from commensuration might be difficult to de- 
tect by diffraction experiments. Similar observation was 
reported for Co-doped BaFe2As2. 6 ' We cannot rule out, 
however, interlayer modulation (1, 0, h + S)o with larger 
S. 

We note that an IC-C crossover has been observed in 
SrFe2As2 under high pressure near 5.4 GPa, where AF 
and SC states coexist, most likely forming nano-scale hy- 
brid structure. 15 ' The 75 As NMR spectra from AF re- 
gion show commensurate peaks at low-T but broad dis- 
tribution at high-T with a crossover near 18 K. Thus 
the incommensurate stripe order may not be exceptional 
but common in iron-based materials, in particular, in the 
vicinity of superconducting phases when Taf is depressed 
below 50 K. 

We next discuss the results in the paramagnetic state. 
The Knight shifts K (= H cS /H) at the 75 As and 23 Na 
nuclei are plotted against temperature in Fig. 4 after 
correcting for the second order quadrupole effects. The 
Knight shifts at the Na sites are small, of the order of the 
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dipolar field from Fe moments. By plotting K at the As 
sites against the susceptibility (the inset of fig. 4), we ob- 
tain the diagonal elements of the hyperfine coupling ten- 
sor: 75 A c = 2.61 ± 0.02 T/fiB, 75 A ab = 3.87 ± 0.07 T//z B . 
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Fig. 4. (Color online) Knight shifts in the paramagnetic state. 
The inset shows the K versus \ plots for the As sites (low-T 
data below 70 K have been omitted since M/H showed upturns 
caused by impurities). 



The lower panel of Fig. 5 shows the temperature 
dependence of the relaxation rate divided by T. The 
Korringa law, (TiT) -1 ~ const., is observed at the 
lowest temperatures both in NaFeAs and in 122 com- 
pounds, 13 ' 14 ) indicating that a part of the Fermi surface 
survives in the AF state even in the presence of a SDW 
gap. The divergence of (T^T)" 1 in NaFeAs near Taf is 
compatible with a second order transition. Then the nu- 
clear relaxation should be dominated by the stripe AF 
fluctuations near the wave vector Q = (1,0, 4) o- The 
anisotropic ratio of the relaxation rate (T-j 110 ') -1 / (T-f) -1 
provides a measure of anisotropy of the AF fluctuations 

14) 

m spin space, ' 

[110k_ 



(7T 



2\S a (Q,uj IC 



\S c (Q,oj res )\ 



(Tf)-i 2|^(Q, Wros )P ' (5) 

where \S l (Q, cj ros )| 2 is the power spectrum of the spin 
fluctuation along the i-direction at Q and the NMR fre- 
quency w ros . The upper panel of Fig. 5 shows diver- 
gent behavior of the ratio (T 1 [11 ° 1 )- 1 /(T 1 C )- 1 of 75 As- 
NMR, which translate into the development of strong 
anisotropy of the AF fluctuations \S a (Q, oj rcs ) | ^ 
\S c (Q,v Tes )\, In BaFe2As2 the divergence is interrupted 
by the first order transition. 

In summary, we have investigated the structural and 
magnetic transitions and AF spin structure in NaFeAs 
using 75 As- and 23 Na-NMR. The NMR relaxation data 
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Fig. 5. (Color online) The lower panel shows the spin-lattice re- 
laxation rate divided by temperature (TiT) - 1 in NaFeAs, and 
BaFeAs 13 ) for two directions of the field. Upper panel shows the 
anisotropy Tf/Tf 10] of the 75 As-NMR. 



prove the second-order nature of the AF transition and 
the strong anisotropy of the AF fluctuations in spin 
space. The most remarkable feature of the AF state is 
the temperature-driven crossover from the high-T incom- 
mensurate structure to the low-T commensurate struc- 
ture. This must be an intrinsic feature of a clean and sto- 
ichiometric system, as demonstrated by the very sharp 
NMR spectra and phase transitions of our crystal. 
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